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Many-body eigenfunctions of the total spin operator can be constructed from the spin and spatial 
wavefunctions with non-trivial permutation symmetries. Spin-dependent interactions can lead to 
relaxation of the spin eigenstates to the thermal eqnilibrium. A mechanism that stabilizes the many- 
body entangled states is proposed here. Surprisingly, in spite coupling with the chaotic motion of the 
spatial degrees of freedom, the spin relaxation can be suppressed by destructive quantum interference 
due to spherical vector and tensor terms of the spin-dependent interactions. Tuning the scattering 
lengths by the method of Feshbach resonances, readily available in cold atomic labs, can enhance 
the relaxation timescales by several orders of magnitude. 
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INTRODUCTION 

A non-degenerate gas of interacting particles far from 
critical points is generally regarded as one of the most 
pronounced representatives of chaotic systems. Accord¬ 
ing to the eigenstate thermalization hypothesis 0,i , ex¬ 
pectation values of observables in the gas eigenstates co¬ 
incide with microcanonical expectation values. The ex¬ 
pectation values relax to the equilibrium after about 3 
collisions, as demonstrated by numerical simulations and 
experiments 

Gases of spinor particles are attracting increasing at¬ 
tention starting from the first experimental |a|7| and the¬ 
oretical ii works (see book [l^ , reviews]o. and 
references therein). Such gases can be described in two 
ways. In the first, conventional, description, each particle 
acquires an additional degree of freedom — the spin pro¬ 
jection Sz, which can have values for spin-i particles. 
It can be either the projection of a real, physical, angu¬ 
lar momentum, or it can be attributed to internal states 
of particles (e.g. hyperfine states of atoms). In the last 
case, the particles can be either bosons or fermions, with 
no relation to their spins. The sum of the particle spin 
projections, the total spin projection, is conserved in the 
absence of spin-changing collisions, being related to oc¬ 
cupations of the spin states. Then the gas is a mixture of 
the gases of particles in the given spin state, which relax 
to the thermal equilibrium with the same temperature. 

Another description of spinor gases is based on collec¬ 
tive spin and spatial wavefunctions. It is a generalization 
of the well-known representation of a two-electron wave- 
function as a product of permutation-symmetric spa¬ 
tial and antisymmetric spin wavefunctions for the sin¬ 
glet state or antisymmetric spatial and symmetric spin 
ones for the triplet state. The singlet and triplet states 
have different energies due to the coulombic interaction 
between electrons. 

The symmetric and antisymmetric functions are exam¬ 
ples of irreducible representations of the symmetric group 
p^l - [l^ . Spatial and spin wavefunctions of V-body sys¬ 
tems with N > 2 can belong to multidimensional, non- 
Abelian, irreducible representations, when permutations 


transform a function to a superposition of the represen¬ 
tation functions. In the case of spin- i particles, the 
representations are associated with the total spin S. The 
total wavefunctions with the correct bosonic or fermionic 
permutation symmetry are expressed as a sum of prod¬ 
ucts of the spin and spatial functions [l3, [13 ■ The only 
one-dimensional representations, the symmetric and an¬ 
tisymmetric functions, are associated with S = N/2 and 
S = 0. Spin-independent interactions between parti¬ 
cles split energies of states with different S, as shown 
by Heitler [I^. The states with well-defined total spin 
are used in quantum chemistry (see Bill) and were ap¬ 
plied to spinor gases [IMl- Many-body entanglement 
of such states can by employed for quantum computing 
. Another example of states with defined spins is the 
Dicke collective state of of two-level atoms or molecules 
coupled by a single mode of the electromagnetic field [s^ 
(see also the recent work [s^ and the references therein). 

In the case of spin-independent interactions, the total 
spin is conserved, the gas can be created in a state with 
given S, and will not relax to the thermal equilibrium, 
which corresponds to a mixture of states with different S. 
The present work analyses the relaxation of such states 
due to spin-dependent interactions between particles. It 
demonstrates that, in spite of coupling to chaotic spa¬ 
tial motion, the spin-relaxation can be suppressed due 
to quantum interference tuned by a Feshbach resonance. 
The relaxation time-scales can be also enhanced in non¬ 
equilibrium ways using time-dependent perturba¬ 

tions. Relaxation of the Dicke states can be suppressed 
due to interaction with cavity modes [s^- 

The paper has the following structure. Section |T] 
describes spin-dependent interactions and permutation- 
symmetric wavefunctions. The Berry’s conjecture 
and eigenstate thermalization Q methods for description 
of the chaos in the spatial degrees of freedom are gener¬ 
alized in Sec. m to the states with well-defined spins. 
In Sec. m these methods are used for calculation of 
non-diagonal matrix elements and relaxation rates, in a 
combination with symmetric group methods and the sum 
rules B,®. 
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I. THE HAMILTONIAN, WAVEFUNCTIONS, 
AND PERMUTATION SYMMETRY 


A general interaction, which does not change the spin 
projection, is a sum of interactions of particles in each 
combination of the two spin states, f and j,, 

Kpin = + al^aVn- (1) 

Whenever the thermal wavelength for the temperature T 

At = [27rftV(wT)]i/2 (3) 

substantially exceeds the interaction range, the interac¬ 
tions can be approximated by the zero-range ones 

% = E t (j))i t (j'))(t (j)i(t (/)i 


vu = E ^ ^ 0’')i 


vn = E t o-))i i (j'))(t (j)ia (/)i 

j¥=j' 

[note the double-counting the particle pairs in Vp and 
Vii, which is compensated by the factors | in Eq. O]- 
The particle coordinates Vj are vectors in H-dimensional 
space (U = 2 or 3). In the two-dimensional (2D) case, 
the motion in the third (axial) dimension is confined by a 
harmonic potential with the frequency Wconf and the two- 
dimensional gas can be formed at sufficiently low tem¬ 
perature T < fiwconf- In certain situations, the two- and 
three-dimensional d-functions should be renormalized in 
order to eliminate divergences. The interaction strengths 


93d — aTTH, -, g2d — 


WWconf ^crcr' 

, (^) 


where m is the boson’s mass, are proportional to the 
elastic scattering lengths aaa-'- 

The interactions Uq-, V^, Vfi, and, therefore, I4pin in 
Eq. ([T]) can be expanded in terms of irreducible spherical 
tensors [ 2 ^ 


Kspin — 2 . 


9DdV+{g]^,-g^,)Vo 
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The spherical scalar interaction 

IgodV = ]-gDd E 




( 4 ) 


spin-dependent parts of the interaction vanish and the 
Hamiltonian of N indistinguishable spin-4 bosons has 
the form 

H = Hq + -godV 

where 

3 

is the kinetic energy and Pj are the momentum operators. 

Since H is invariant over independent permutations of 
the particle spins and coordinates and commutes with 
the operators of the total spin S and i ts proj ection Sz, 
the eigenfunctions can be expressed as [IJ, [II, [^ 

= (5) 


(S) . '-^(S) 

where the spatial and spin wavefunctions form 
bases of irreducible representations of the symmetric 
group of permutations of N symbols [I3l - [l5l |. The 
representations are associated with the two-row Young 
diagrams A = [N/2 + S, N/2 — S] and have dimensions 

Ai!(25 + 1) 

~ {N/2 + S + iy.{N/2 - sy: 

The basis functions within the representations are labeled 
by the standard Young tableaux t of the shape A. A per¬ 
mutation V of the particles transforms each function to 
a linear combination of functions in the same represen¬ 
tation, 


t' t' 

( 6 ) 

Here D[^j (V) are the Young orthogonal matrices. Their 
properties [l3l - ll5] provide the correct bosonic transfor¬ 
mation for the total wavefunction ([5]). 

The explicit form of the spin wavefunctions [s^ is not 
used here. Their orthonormality 




Ss'sSt’tSs’^s^ 


leads to the Schrodinger equation for the spatial wave- 
functions 

(7) 

(all wavefunctions within an irreducible representation 
are energy-degenerate, according to the Wigner theo¬ 
rem). 


H. QUANTUM-CHAOTIC WAVEFUNCTIONS 
WITH DEFINED TOTAL SPINS 


with the interaction strength g^d = {god + 9Dd + 
provides the spin-independent interaction between par¬ 
ticles. If all scattering lengths have the same value, the 


Consider N spin-4 bosons in a periodic box with in¬ 
commensurable dimensions. The box can be either three- 
dimensional (3D) of the volume or 2D of the square 
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with tight confinement axial dimension. A tight con¬ 
finement in two directions can lead to a homogeneous 
one-dimensional gas (see [i^), which is an integrable 
system and is not considered here. In contrast, homo¬ 
geneous 2D and 3D gases can demonstrate chaotic be¬ 
havior at the sufficiently high energy-density of states, 
as assumed in the first consideration of eigenstate ther- 
malization by Deutsch [l|. According to the quantita¬ 
tive criteria |4ll44^ . based on analyses of delocalization 
in the Fock space, many-body systems become chaotic 
when the interaction matrix elements exceed the energy 
spacing of directly-coupled many-body states. In the 3D 
case, the matrix elements are g^d/L^- The two-body in¬ 
teractions couple states with different momenta of any 
of N{N — l)/2 pairs of particles and conserved center- 
of-mass momentum. Then the energy density of coupled 
states will be N{N — l)/2 times the energy density of 
relative-motion states L^'rvP'I'^T^I'^j This leads 
to the criterion of chaos 

_ arr + “ii + “tt ^ v^At 
as =- 3 -> ^2 ^ (8) 

where the thermal wavelength Ay is given by Eq. ©■ 
This criterion do not contradict to the sufficient condi¬ 
tion 0 for the hard-sphere gas, based on calculations 
for chaotic billiards, that chaos appears when the parti¬ 
cle radius exceeds the thermal wavelength. Compared to 
this condition, the present criterion reduces the thresh¬ 
old temperature by the factor . In the 2D case, the 
energy density of relative-motion states is mLi^/{A ttIt?) 
and the matrix elements are g 2 dlL'^- Then the condition 
of chaos will be 

as > 2 \/27rft/ (mwconf(9) 

where the square root is, up to a factor, the confinement 
range. 

Under the conditions of chaos, the spatial wavefunction 
can be represented according to the Berry conjecture [s^ . 
In the Srednicki form , the non-normalized solution of 
Eq. ([7]), labeled by the index n, with the well-defined 
total spin S' is a superposition of plane waves 


cx ^ A(f)(<,{p})J({p}2 - 2mEi^^) 

{pl 


X exp(i 

3 


( 10 ) 


with the momenta Pj in the periodic box with incom¬ 
mensurable dimensions. Due to discrete spectrum of Pj, 
the states with approximately fixed energies En ' are se¬ 
lected by the function 

Six) = e{A-\x\)/i2A), 


where 0(x) is the Heaviside step function. In the final 
calculations, when the summation over {p} is replaced 
by integration, (5(a;) is replaced by the Dirac (5-function. 


Then the total wavefunction m can be represented as 
(see Appendix lAl) 

^ E {p})5({p}2-2mE(^))^'(fJ,j5 

{p} 

( 11 ) 

It is a superposition of symmetrized plane waves — wave- 
functions of non-interacting particles (see Eq. 

(ESI) and [l^.ffalrijll. Given S and Sz, these wavefunc- 
tions are labeled by the Young tableau r and the set of 
particle momenta {p} = {pi,... ,PAr} ({p}^ = P])- 

The summation over the simplex Pi < P 2 < ■ • • < Pat 
is denoted as , where p < p' if or = p(, 

and py < p'y, or p^, = p^, py = p'y, and Pz < p'z- Such 
summation, neglecting multiple occupations of the mo¬ 
mentum states, is applicable to non-degenerate gases, 
when the difference between Bose-Einstein, Fermi-Dirac, 
and Boltzmann distributions is negligibly small. The 
normalization factor (see Appendix ICl) provides 

(<sj<sl) = l- 

According to the Berry’s conjecture [1, the co¬ 
efficients A(f^(r, {p}) are treated as Gaussian random 
variables with a two-point correlation function (see Ap¬ 
pendix IbJ , generalizing the one of 0 to the states with 
well-defined spins. 


„ V „ V ^({p/}2 _ {p}2) 

( 12 ) 

Here, as in 0, ()gg denotes average over a fictitious 
“eigenstate ensemble”, which describes properties of a 
typical eigenfunction. The Kronecker s;rabols appear 
here, as well as in the correlation function (2|, since differ¬ 
ent S and n correspond to different (independent) eigen¬ 
functions and different {p} in the same simplex corre¬ 
spond to different (independent) plane waves. In ad¬ 
dition, Eq. (HU contains the Kronecker symbol of the 
Young tableaux r and r', as proved in Appendix [BJ 


III. DECAY RATES 


The rate of transitions from the state with the spin 
S to the S' one is estimated by the Weisskopf-Wigner 
width (see 0 ) 


(13") 

where the density of states /dE is evaluated in Ap¬ 
pendix o Eor a typical wavefunction (HU, the squared 
modulus of the matrix element can be estimated by the 






spml'k^g. 




dE 
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eigenstate-ensemble average 

^({P}^ - ‘2rnEi^^)d{{p'y - 2mEl^ 

{p} {p'l 

X K^r'-rLl.S' l^spinl^rlpr.?,)!^- (1^) 

rr' 

Here S' ^ S and the product of four coefficients 
4ln^^(r, {p}) leads to a four-point correlation function 
(|B5l) . which is represented by Eq. (|B4I) in terms of two- 
point ones m for the Gaussian ensemble as in Q . This 
relation of matrix elements between wavefunctions of in¬ 
teracting and non-interacting particles is obtained since 
the correlation function contains Sr'r- The sum of 
squared moduli of the matrix elements between the wave- 
functions of non-interactin g pa rticles in Eq. (HI is cal¬ 
culated with the sum rules |28l |. 

The expansion (jT]) of the interactions 14pin contains ir¬ 
reducible spherical scalar, vector, and tensor. According 
to the Wigner-Eckart theorem, the scalar interaction V 
conserves the total spin, while the spins S and S' can 
be coupled by the vector component Vq if [S' — 5"| < 1 
and by the rank 2 tensor component Vg ' if [S' — S"| < 2. 
This leads to quantum interference of the vector and ten¬ 
sor contributions to the transitions between states with 
spins S and S' = Si 1. Equation (ITT)) and the sum rules 
lead (see Appendix IdI) to the transition rates 


dAS-i) _ {S^-S^,){N + 2S + 2) 


S{2S + 1)N 


S^AN + 2) 

■ 52-1 


+ ai{N -2)+4:a+a-S^ 


r_Dd 

(15a) 


dAS+i) ^ [iS+l)^-S^J{N-2S) 
(5'-bl)(2S'-tl)A^ 


■ S{S + 2) 


+ a^_{N -2)+Aa+a_S^ 


r_Dd 

(15b) 


p(5.S-2) ^ [(^-l) 2 - 52 ](^ 2 _g 2 ) 

2S{S-1){2S-1){2S+1)N 

x{N + 2S){N + 2S + 2)alTDd (15c) 

p(5.S+2) ^ [(5+l)2-,g2][(g + 2)2-g2] 

2{S + 1){S + 2){2S + 1){2S + 3)N 

X {N-2S){N-2S-2)alTDd (15d) 

where a+ = (a^i Ea^- 2a^^)/as, «_ = (o^i - a^)/as, 
as is given by Eq. m, and the interference terms are 
proportional to a+a_. 



S 

FIG. 1. The total decay rates for the state with the total spin 
S and its projection Sz for two values of S — Sz- The solid 
black lines correspond to the background scattering lengths. 
The dashed blue and dot-dashed red lines show the decay 
rates minimized by tuning of an and respectively. 


The rates are proportional to the factors (see Appendix 

IQ 


Tsd = 2^J ttT ! masn^d, ^2d = ^a|wconfn2d (16) 

where nud is the H-dimensional gas density and the 
temperature is related to the eigenstate energy as T = 
2Ei^^/{3N) 0. In the 3D case (D = 3), Ta^, up to a 
numerical factor, is the frequency of elastic collisions per 
particle in the gas. In the 2D case, T 2 d is proportional 
to the rate of collisions per particle too, since the proba¬ 
bility of collision during one axial oscillation is 87ragn2d, 
the oscillation frequency is 27ra;conf, and the oscillation 
velocity substantially exceeds the 2D motion one in the 
2D regime {T < Swconf)- The present derivation is valid 
whenever T substantially exceeds the degeneracy tem¬ 
perature Tdeg = 2TTh^nyj^/m [l^ . 

For the two states of ®'^Rb atoms , generally used 

in experiments, | 4-) = IT’ = I,to/ = I) and | 

) = |F = 2,m/ = —I), the scattering lengths are [451 
nu ~ 100.4as, aq ~ 95.5os, and a't- 4 , ~ OS.Oob, where 
as is the Bohr radius. For 713 ^ = 10^^cm“^ we have 
Tdeg ~ 0.04/rK. The zero-range approximation ([3]) is 
applicable whenever T <C 10"^/iK and, according to the 
criterion m, the system of A^ = 10 ^ atoms becomes 
chaotic at T > 8 X 10“^^/iK. Then for T = 1/iK we 
have Tad ~ 0.9s“^. 

The total decay rate Tg^^ = Ss'/s ^ presented 
in Fig. [H Even for Af = 10^ atoms it does not exceed 
~ GOr^d- The decay is suppressed at large values of S 
and Sz- 

The Feshbach resonance tuning of the elastic scatter¬ 
ing length [ 4 ^ cannot eliminate the decay, as one tuning 
parameter — the magnetic field — cannot make van¬ 
ish two combination of the scattering lengths — and 
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FIG. 2. The total decay rates for the state with the total spin 
S and its projection Sz minimized by 0 - 1 - 4 , tuning in the vicinity 
of the Feshbach resonance at Bo « 9.13G in ®^Rb. The solid 
black, dashed blue, and dot-dashed red lines correspond to 
S — Sz ~ 100, 20, and 5, respectively. The plots of optimal 
detuning B — Bo for different S — Sz almost coincide. 



S 

FIG. 3. The total decay rates for the state oi N = 10"^ atoms 
of ®^Rb with the total spin S and its projection Sz minimized 
by 044 tuning in the vicinity of the Feshbach resonance at 
Bo « 1007.4G. The solid black, dashed blue, and dot-dashed 
red lines correspond to S — 5'z = 100, 20, and 5, respectively. 
The optimal detuning B — Bo for different S — Sz is presented 
by three upper plots. 


a_. However, the Feshbach tuning can minimize the 
decay due to the destructive interference of the contri¬ 
butions of the spherical vector and tensor interactions, 
mentioned above. These contributions are proportional 
to a- and a+, respectively, in Eq. (ITSl) . In the case of 
®’^Rb, the resonance in fi collision at Bq ^ 9.13G (with 
the width A Ri 15mG) is well separated from the res¬ 
onances in collisions at Bq > 400G (the widest one 
at Bq R! 1007.4G has A r 0.21G) [i^. No resonances 
are known in ff collisions. Then either a-t- 4 , or 044 can 
be tuned without changing other scattering lengths. The 


resulting decay rates are presented in Fig. [T] Tuning of 
044 can reduce the decay rate at small S, while tuning 
of a -44 can lead to the reduction by orders of magnitude 
at large S and Sz- The magnetic field detunings, mini¬ 
mizing the decay rates, and the minimal rates are plot¬ 
ted in Figs. [Hand [21 Even the minimal detunings 0.5G 
and 4G, respectively, substantially exceeds the resonance 
widths. Then the closed-channel effects can be neglected, 
although the resonances are closed-channel dominated 
( 4 ^ . The resonance-enhancement of three-body losses 
(see [13) can be neglected at such detunings as well. 

Approximate expressions can be obtained in the re¬ 
gions of maximal suppression, Sz ~ iS. Here the a -44 
tuning minimizes the decay rate when a+ R To?- and 
the minimal rate is approximated for A ^ 1 by 

- - (A - 2S)N{4S + 5)] |. 

Then for A/2 — S = const <C A, S' — |Sz| = const, and 
fixed density the decay rate is scaled as I/A. The desired 
scattering lengths, a -44 ~ ®TT ~ ®T-1 ~ 

Sz ~ —S , are obtained for ®^Rb at B — Bq fv 0.6G or 
— 0.6G, respectively. The state with S = Sz = A/2 for 
even A (or with S = Sz = (A — I)/2 for odd A ) cannot 
decay to states with S' < S since S' cannot be less than 
Sz, nor to states with S' > S since S' cannot exceed A/2. 

When the decay rates are substantially suppressed, 
the lifetime of the states with well-defined spins is re¬ 
stricted by the loss processes in the cold gas, such as 
spin-changing (dipolar relaxation) and three-body colli¬ 
sions, leading to high-energy atoms escaping the trap. 
The spin-changing collisions can also change many-body 
spin and its projection, leading to additional decay. How¬ 
ever, this decay, having the rate comparable to the gas 
loss rate, does not lead to additional restriction of the life¬ 
time of the states with well-defined spins. In real physical 
situations, dipolar relaxation becomes substantial only 
for Cs and atoms with high magnetic momenta [HI- 

CONCLUSIONS 

Spatially-chaotic many-body eigenstates of the total 
spin operator can be described, according to the Berry 
conjecture [s^, within the Srednicki approach 0 [see Eq. 
mi. This description, in a combination with the sum 
rules [13, allows us to evaluate the matrix elements of 
spin-dependent two-body interactions, leading to transi¬ 
tions between states with different total spins. The tran¬ 
sition rates, calculated within the Weisskopf-Wigner ap¬ 
proach m, are proportional to the elastic collision rate 
per particle (HU. The decay rates can be suppressed 
due to destructive interference of the contributions the 
spherical vector and tensor terms in the spin-dependent 
interaction V/pin. The interference terms in Eq. m can 
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be controlled by Feshbach resonances as they are propor¬ 
tional to a_|_a_. Another manifestation of quantum inter¬ 
ference is the effect of dynamical localization. It can slow 
down relaxation due to periodic driving while 

the present mechanism pertains to time-independent sys¬ 
tems. The long-lived entangled states can find applica¬ 
tions in quantum computation and metrology. 


Appendix A: Wavefunctions with defined total spins 
of interacting and non-interacting particles 


A permutation of coordinates in the wavefunction m 
can be represented as 

cx ^Ai^)(t,{p})5({p}2 - 2mEify) 

{p} 

j 

cx ^ Af)(t,P{p})5({p}2 - 2mF;f)) 

{P} 


X exp(l 

j 


are spatial wavefunctions of non-interacting particles [l^ . 
[TsI . . They satisfy relations ([HI) , forming a basis of the 
irreducible representation associated with A. The Young 
tableau r labels different representations for the same 
S and {p}. The total wavefunction of non-interacting 
particles 

(A3) 

t 

is expressed similarly to [see Eq. ([5])]. Then Eq. 

(IA2I) leads to Eq. (HU. The wavefunctions of non¬ 
interacting particles form an orthonormal basis 

= Ss'sSs'^sJr'rS^^'n^} (A4) 


and satisfy the Schrddinger equation 


r{p}s, - ^np}s.- 


Appendix B: Two- and four-point correlation 
functions of the coefficients {p}) 


where 'P{p\ = {p-pi,..., Pvn}- Then Eq. ([5]) leads to 
A(^) (t, V{p}) = ^ 4^' (P) A(^) (r, {p}) (Al) 

r 


According to the Berry conjecture the coefficients 


(t, {p}) in Eqs. (ITUl) and (fTTl) can be treated as Gaus¬ 
sian random variables with a two-point correlation func¬ 
tion 0 


( S) 

Let us represent the wavefunction in the form 
which explicitly shows its permutation properties. This 
can be done by summation in Eq. m over the simplex 
Pi < P 2 < • • • < Pjv- Momentum sets in other simplices 
are given by 7^{p}. Neglecting contributions of the sets 
{p} which contain equal momenta Pj = Pji, one gets 

cx ^^'A(f)(t,P{p})5({p}2 - 2mE;(^)) 

-P {p} 

X exp(i ^ p'pjTjjh) 
j 

« E' E ^ (a {p})^({p}" - ^ (V) 

{p} P 

X exp(i^pjrp-ij/ft) 


_ Ss'S^n'nS{p'}{p}f{t,t') 

EE " 5 ({p '}2 _ {p} 2 ) 

(Bl) 

where (t, {p}) = {— p}). The Kronecker 

symbols appear here since different S and n correspond to 
different (independent) eigenfunctions and different {p} 
within the given simplex correspond to different (inde¬ 
pendent) plane waves. (The correlation functions with 
{p} and {p'} in different simplices do not appear within 
the present paper.) By now nothing can be told on 
since and are components of the same 
eigenfunction, related by Eq. (|6]). The factors are 

determined below. 

Equation m leads to the following equality for an 
arbitrary permutation V 


(Aif)*(f',{p'})Ai^)(t,{p}) 


J 

« Y.'Y.^nHr,{Pm{Pr-2rnEif^) 

{p} 


-( 5 ) A(^)(t,{p}) =^ZAW(P)A(«)(r,P{p}). 

[hJ 

(A2) It relates coefficients Ai‘^^(r, {p}) in different simplices 
and leads to 


Here Eq. m and the identity for the orthogonal repre¬ 


sentation matrices d]^^(V ^) = d\^\v) are used and 

^MpI = (^) E i'P) exp(f ^ Pjrvj/K) 

\ ■ / -p j 


(A(«)*(f',{p})A(«)(t,{p})) 

\ / hjhj 

= {Ai^^*{r',V{p})Ai^\r,np})) 

\ / Ilill/ 
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Applying Eq. dm to all correlation functions, we get 


one gets the equation 


fit, t') = ^ d]^1, W iV)fir, r') (B2) 

rr' 


for any V. Averaging the right hand side of this equation 
over all V and using the orthogonality relation [IJ, Il5| 




Its solution is fit,t') = constJtt'. Equation (fT^ is ob¬ 
tained if const = 1. Other choice of the constant factor 
can only change the normalization factor in Eq. dill) . The 
function fit,t') = 5tt' satisfies Eq. (jB2p . as can be easily 
proven using the identity 

^ D^^liV)D^^} iV) = iVV-^) = 5,., (B3) 

r 


= —Su'Srr'S^y, 


for the Young orthogonal matrices (see 0^3). 

The relation between four-point and two-point corre¬ 
lation functions 


EE 


\ ' ibrj \ ' \ 

+(A^;')(t",{p"})Air)(t',{p'}))^^ ^Af,:\t"',{p"'})A^^\t,{p\))^^ (b4) 

is a straightforward generalization of the similar relation in Q. For the product of four coefficients {p}) in 

Eq. m it leads to 


EE 


' EE 


\n \ , ll' JJ n \ ^ n \ , ll' J J n \ > ^VeE d('In"I2 - Ir)I2U('/n'"I2 - Ir)'I2'> ^ ^ 


r'Sr 


EE d({p"}2 - {p}2),5({p'"}2 - {p'}2) 


since S' S. 


Appendix C: Normalization factor and the density 
of states 

Orthonormality of the non-interacting particle wave- 
functions (lAH) leads, using the two-point correlation 
functions to the following overlap of the wavefunc- 
tions of interacting particles (EH) 

fsSs'^sJs'sSn'u 

X (Cl) 

{p} 

The summation over {p} in the simplex can be approx¬ 
imated by integration over whole momentum space with 


the replacement of S by the Dirac i5-function, 

{p} ^ ^ 

X J d^^p6i{p}^-2mEif^) (C2) 

The integral 

/d'=V({pr-£) = h|gT ,C3) 

is calculated in [3 . Then the normalization factor is given 
by 

2Nl{2TT)D^/^r{DN/2)h^^' ^ ’ 

The density of states of interacting particles can be ap¬ 
proximated by the one of non-interacting ones. Number 








of such states with the total spin S below the energy E 
is 

n^^HE) = - {P}V(2m))- 

{P} 

Then the density of states is estimated as 
(E) ^ ^ iE +A)- n^^'> (E - A) 

= fsJ2'~^iE - {p}V(2m)) = 2m )) ”' (C5) 

{P} 

[cf. Eq. ([CT])]. 

Appendix D: Matrix elements and decay rates 


elements (p'p|(5|pp') from their average values [1^ and 
vanish in the present case since 

(p'p',|d|pjPjv) = A“'^^p'.+p',p,+p^., 

and the matrix elements (p^p|d|pp') = L~^ are con¬ 
stant. The sets {p} and {p'} different by single momen¬ 
tum are not coupled due to the momentum conservation. 

The form (ID2p of the sum implies that the unchanged 
momenta are in the same positions in the sets {p} and 
{P'}- However, arbitrary permutations V and V' of the 
momentum sets do not change the sum, 

r,r' 

r,r' 

(D3) 


If the scattering lengths and interaction strengths 
are spin-dependent, states with different total spins 
become coupled by the spherical vector Vq and tensor 
interactions in Eq. ®. Matrix elements of the 
spherical vector and tensor are related to ones for the 
maximal allowed spin projections by the Wigner-Eckart 
theorem 


It/ \_ n/ 

\_ y/flyl'S'O iT/ 1^1 


(Dl) 


and do not vanish for \S — S"| < 1 and \S — 5"| < 2, 
respectively. Here the ratios of the 3j Wigner symbols 
^S^’o tabulated in [23, [l^- Matrix elements with 

S' > S are calculated using Hermitian conjugate in Eq. 
(ID1|) . taking into account that (Vo)’l' = Vq, = 

(see [ 2 ^). 

Sums of the products of the matrix elements of the 
spherical tensors between the wavefunctions of non¬ 
interacting particles (IA3I) can be represented for S' < S 

as [111 




r'{p'}S 






Vs' 

N(N -1) 


E l(PjPTl'5|PjPi')l^ 

3 < 3 ' 


3 '¥^ 3 " ¥^3 


Here the universal factors 14 ] are independent 

of the occupied spatial modes, while the sum of matrix 
elements is independent of the total spin and its projec¬ 
tion. The sum rule (|D2p is valid if the sets {p} and 
{P'} are different by two momenta. Similar sum for 
{p'} = {p} is proportional to deviations of the matrix 


This equality is provided by the transformation of the 
non-interacting particle wavefunctions on permutation of 
the quantum numbers (see 0 ) 




(S) 

rP{p}S 




(S) 

r'{p}S 


r' 


and the identity (IB3D . Then the sum rule (ID2I) can be 
rewritten as 


E(^ 


(S') 


,{S) 


r^pqS^ l^l^r{p}S 




r{p}S 




= Y^^’^'>[Va,Vb] 


fs' 

N{N- 1) 


E E l(PpjPpj'l'^IPiPT)l^ 

V j<j' 


^ n '^p^j/'Pj"' 

3 '¥ 3"¥3 


The factor 2 in Eq. (ID2]) is absent here since the Kro- 
necker symbols are satisfied by two permutations, V and 
Wjj>, when the momenta Pj and Pj' are changed. 

In the matrix element (II4p between the wavefunctions 
of interacting particles, we can replace the sum over {p} 
in the simplex pi < p 2 < • • • < Pat by average of sums 
over all simplices, Pvi < PP2 < ■ ■ ■ < Pvn, 

EA{p}) = Ay:y:V(P{p}). 

{p} ■ -p Pfp} 


Since the sums of the matrix elements (lD4l) and sums 
of squared momenta {p}^ are invariant over momentum 
permutations, Eq. (IT4)) takes the form 



X E"E"^({P}' - 2mE(^))d({p'}2 - 2mBif)) 

{p} fp'l 

r,r' 
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where X]{p} denotes summation over the sets {p} which 
do not contain equal momenta. 

The Weisskopf-Wigner decay rates m for transitions 

from the S'-multiplet to the S' one, \ calculated 
with the matrix elements (ID5I) of the operator ([4]), rela¬ 
tions m, and sum rules (ID4I) attain the form 


qs,s-i) 




4] 




fs-i 

J Nfs 


(D6a) 


.(S.S+l) 


^ V(^+1’2)[4), 1/4] 


(D6b) 


(HU), dns]), and (lD4l) as 

TTfsNgl, ( MPMp V dn^^'HEP) 

2hL^D iV! j dE 

^ - 2to4^4({p'}^ - 2m4f 

{p}#{p'} 

x^i^EE 4i+Ppj'Pi+P3' n 

^ ' j<j' V 

The Kronecker symbols here fix values of all p', except for 
p' — p',. The relation (IC5I) between the density of states 
and the normalization factor and replacement summation 
by interaction lead to 


Tud = 


2 / T \ D(N+1) 

'^’mgpd ( L \ 

2 ftL 2 Z>^! \2T:h) 




c J p6{{p}^—2mE^^) J d^p'6 ^2p'^ — 


(Pi - P2y 


Calculating the integrals over pa,... ,pjv and pi -f p 2 
with Eq. (IC3|) and using Eq. (jC4l) for the normalization 
factor, we get 


Tcd 


_ jVr(j^iV/ 2 )mg 4 _ 

2^+27r3^^/2-ir(i:i(iV - l)/2);i^^+iL^(mE4)^^/^”^ 



r('S,S- 2 ) 




-2,V-lh^^Dd 


'Nfs' 
(D6c) 


Since f d?p'5 (^p'^ — p^^ = e and J d^p'S (^p'^ — P^^ = 
2ttp, integration over p leads to Eq. ®. 

Substitution the factors 4^0^ ’"^^ ^’^d T('^’2)[t4,14] 
from [13, [13 into Eqs. (ID6|) leads to Eq. m- 


-.(3,3+2) 


2 / 

5 (a.V 


(3+2,3. 
3^0 


'^^yY(3+2-2)[v(^lv^^)]lT 


Dd 


'TV 
(D6d) 

Here the factor Tpd is independent of S, S', and Sz, as 
will be explicitly seen below. It is obtained from Eqs. 
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